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Abstract 

All possible hermitian texture 6 zero Fritzsch like as well as non Fritzsch like 
lepton mass matrices (144 combinations in all) have been investigated for both Ma- 
jorana and Dirac neutrinos for their compatibility with the current neutrino oscilla- 
tion data, keeping in mind the hierarchy of neutrino masses. All the combinations 
considered here are ruled out by the existing data in the case of inverted hierar- 
chy and degenerate scenario of Majorana neutrino masses. For Majorana neutrinos 
with normal hierarchy, we find that out of 144 combinations, only 16 possibilities 
can accommodate the experimental data. Assuming neutrinos to be Dirac particles, 
normal hierarchy, inverted hierarchy as well as degenerate neutrinos are ruled out 
for all the combinations. 



1 Introduction 

Evidence of neutrino oscillations, as obtained from solar neutrino experiments PQ, at- 
mospheric neutrino experiments [2], reactor based experiments [3] and accelerator based 
experiments jl], has firmly established the existence of massive neutrinos and mixing of 
lepton flavors. In the last few years, there has been considerable progress in the mea- 
surement of neutrino mass square differences and mixing angles, however, the absolute 
neutrino mass scale and the smallest lepton mixing angle 6*13 are still unknown. Further, 
the presently available neutrino oscillation data does not throw any light on the neutrino 
mass hierarchy, which may be normal/inverted hierarchy and may even be degenerate. 
Furthermore, the situation becomes complicated when one realizes that neutrino masses 
are much smaller than charged fermion masses as well as it is yet not clear whether 
neutrinos are Dirac or Majorana particles. In the absence of a convincing fermion flavor 
theory in this regard, several approaches have been adopted, e.g., texture zeros [5], seesaw 
mechanism [6], radiative mechanisms [7], flavor symmetries [8], extra dimensions pJJ, etc.. 
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In this context, texture specific mass matrices have received a good deal of attention in 
the literature, in particular, Fritzsch-like texture specific mass matrices seem to be very 
helpful in understanding the pattern of quark mixing and CP violation [101 HI] . Taking 
clues from the success of these matrices in the context of quarks, several attempts have 
been made to consider texture specific lepton mass matrices [12]- [IS] for explaining the 
pattern of neutrino masses and mixing. In the absence of a sufficient amount of data 
regarding neutrino masses and mixing, it would require a very careful scrutiny of all pos- 
sible textures to find viable structures, which are compatible with data and theoretical 
ideas, so that these are kept in mind while formulating mass matrices at the grand unified 
theory (GUT) scale. 

In the quark sector, both Fritzsch like as well as non Fritzsch like texture six zero mass 
matrices have been completely ruled out [17] . In leptonic sector, most of the analysis have 
been carried out in the flavor basis [H], whereas in non flavor basis, when both charged 
lepton mass matrix M; and neutrino mass matrix M v are three zero type, the number 
of possibilities for texture six zero mass matrices becomes very large. These possibilities 
have been explored in the literature [T5] both for Majorana as well as for Dirac neutrinos, 
however adequate attention has not been given to the cases for inverted hierarchy and 
degenerate scenarios. Also, it is perhaps desirable to note that Dirac neutrinos have not 
yet been ruled out by the experiments [18], therefore it becomes important in the case 
of texture six zero possibilities to carry out a detailed comparison for Dirac like as well 
as Majorana like neutrinos in the case of normal, inverted and degenerate cases. This 
exercise becomes all the more interesting in view of the refinements of data and advocacy 
of quark lepton symmetry [19J. 

The purpose of the present paper is to update the analysis of [15] for all possibilities 
of texture six zero lepton mass matrices as well as to extend this analysis to the case 
of inverted and degenerate neutrino masses. To preserve the parallelism between quarks 
and leptons only those neutrino mass matrices have been considered, which are consistent 
with the requirement of non zero and distinct neutrino masses. Following our analysis in 
the quark sector [17J, the mass matrix for leptons and neutrinos are taken to be hermitian. 
For the sake of completion, we have also investigated the cases corresponding to charged 
leptons being in the flavor basis. It would also be a desirable exercise to calculate several 
phenomenological quantities, such as effective neutrino mass (m ee ) related to neutrinoless 
double beta decay, Jarlskog's rephasing invariant parameter in the leptonic sector Ji and 
the corresponding Dirac like CP violating phase S\ for the viable cases. 

The detailed plan of the paper is as follows. In section 2, we present the essentials of 
the formalism connecting the mass matrix to the neutrino mixing matrix. Inputs used 
in the present analysis have been given in section 3. In section 4, various possibilities of 
texture six zero mass matrices have been given. In sections 5 and 6, for Majorana and 
Dirac neutrinos respectively, the detailed calculations pertaining to normal/inverted and 
degenerate cases have been discussed. Finally, section 7, summarizes our conclusions. 
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2 Construction of PMNS matrix from mass matrices 



To begin with, we present the Fritzsch like hermitian texture six zero lepton mass matrices, 
e.g. ; 

/ A t \ / A v \ 

M { = [At Bi , M vD — I A* v B v , (1) 
\ Bf Q ) \ Bl C v ) 

Mi and M v d respectively corresponding to charged lepton mass matrix and Dirac-like 
neutrino mass matrix. It may be noted that each of the above matrix is texture three 
zero type with = \A^ U ) |e* ai W and B^ = \B^ u )\e ll3l< -^ . For Majorana neutrinos, the 
neutrino mass matrix M v is given by seesaw mechanism |6], for example, 

M v = -M T vD {M R y x M vD , (2) 

where M v d and Mr are respectively, the Dirac neutrino mass matrix and the right-handed 
Majorana neutrino mass matrix. It may be mentioned that for both Majorana as well 
as Dirac neutrinos the texture is imposed only on M V Bi with no such restriction on M v 
for Majorana case. In the absence of any guidelines for the right handed Majorana mass 
matrix Mr as well as to keep the number of parameters under control, it would be desirable 
to keep its structure as simple as possible, therefore we take Mr = turI where I is the 
unity matrix and itlr denotes a very large mass scale. 

To fix the notations and conventions as well as to facilitate the understanding of 
inverted hierarchy case and its relationship to normal hierarchy case, we detail the for- 
malism connecting the mass matrices to the neutrino mixing matrix. To facilitate the 
diagonalization of M^, where k = I, vD, the complex mass matrix M^ can be expressed 
as 

M k = Q k M r k P k (3) 

or 

M r k = Q\M k P\ , (4) 

where Ml is a real symmetric matrix with real eigenvalues and Q k and P k are diagonal 
phase matrices. For the hermitian mass matrix Q k = P^. In general, the real matrix Ml 
is diagonalized by the orthogonal transformation O k , e.g., 

M Ma 9 = T M r 0k j (5) 

which on using equation (j3J) can be rewritten as 

M* ag = O k T Q[M k PlO k . (6) 

To facilitate the construction of diagonalization transformations for different hierarchies, 
we introduce a diagonal phase matrix defined as diag(l, e l7T , 1) for the case of normal 
hierarchy and as diag(l, e t7T , e in ) for the case of inverted hierarchy. Equation ([6]) can now 
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be written as 

£ fe M? a9 = O k T Q\M k PlO k , (7) 
which can also be expressed as 

Mt 9 = iWQlM k PlO k . (8) 
Making use of the fact that 01 = O k it can be further expressed as 

M d k ia9 = (QkO k £ k yM k (PiO k ), (9) 

from which one gets 

M k = Q k O k i k Mt ia9 T k P k . (10) 

The case of leptons is fairly straight forward, whereas in the case of neutrinos, the 
diagonalizing transformation is hierarchy specific as well as requires some fine tuning of 
the phases of the right handed neutrino mass matrix M R . To clarify this point further, 
in analogy with equation (flOl) . we can express M v d as 

M vD = Q vD O vD i vD M^ 9 T vD P vD . (11) 

Substituting the above value of M vD in equation (j5J) one obtains 

M v = -{Q vD O vD t vD M v h * 9 Ol D P vD ) T {M R y\Q vD (12) 

which, on using Pj D = P u d, can further be written as 

M v = -Pu D O uD M^ 9 i uD Ol D Ql D {M R )- l Q uD O vD i uD M^ 9 Ol D P uDl (13) 

wherein, assuming fine tuning, the phase matrices Q^ D and Q v d along with — M R can be 
taken as m R diag(l, 1, 1) as well as using the unitarity of £ U £> and orthogonality of O u d, 
the above equation can be expressed as 



M v = P vD O vD y uD ' Oi D P vD . (14) 
m R 

The lepton mixing matrix or the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix 
[20] U can be obtained from the matrices used for diagonalizing the mass matrices Mi 
and M v and is expressed as 

U=(QiO&y{P vD O vD ). (15) 

Eliminating the phase matrix by redefinition of the charged lepton phases, the above 
equation becomes 

U = 0}Q l P 1/D O uD , (16) 

where QiP u d, without loss of generality, can be taken as diag(e** 1 , e^ 2 , e^ 3 ); <pi, 2 and 
03 being related to the phases of mass matrices and can be treated as free parameters. 
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When Mi and M v both are texture three zero type, it is sufficient to consider only two 
phases, i.e. <pi and <p 2 or <pi and <p 3 or (p 2 and <p 3 . 

3 Inputs used in the present analysis 

Before going into the details of the analysis, we would like to mention some of the essentials 
pertaining to various inputs. The inputs for masses and mixing angles used in the present 
analysis at 3a C.L. are as follows [21J, 

Am| 2 = (6.90 - 8.20) x 10" 5 eV 2 , Am 2 3 = (2.09 - 2.83) x 10~ 3 eV 2 , (17) 

12 = 31.7° - 37.7°, 9 23 = 35.2° - 53.7°, 13 < 12.7°. (18) 

The above data reveals that at present not much is known about the hierarchy of neutrino 
masses as well as about their absolute values. 

The masses and mixing angles, used in the analysis, have been constrained by the 
data given in equations ffTTj) and f|T8|) . For the purpose of calculations, we have taken 
the lightest neutrino mass and the phases <pi and <p 2 as free parameters. The other two 
masses are constrained by Am 2 2 = m 2 2 — m 2 ^ and Am 2 3 = m 2 3 — m 2 2 in the normal 
hierarchy case defined as m Ul < m U2 <C m V3 and also valid for the degenerate case defined 
as m vi m v2 ~ m vi an d by Am% 3 = m^ 2 — vr? V3 in the inverted hierarchy case defined 
as m U3 <C m Ul < m U2 . It may be noted that lightest neutrino mass corresponds to m Ul 
for the normal hierarchy case and to m V3 for the inverted hierarchy case. The explored 
range of lightest neutrino mass is taken to be 0.0001 eV — 1.0 eV as our results remain 
unaffected even if the range is extended further. In the absence of any constraints on the 
phases (pi and (p 2 these have been given full variation from to 2n. 

4 Texture six zero lepton mass matrices 

To begin with, we enumerate the number of possibilities for the texture six zero lepton 
mass matrices. It is easy to check from equation (TjQ) that there are 20 possible patterns of 
texture three zero hermitian mass matrices, which differ from each other with regard to 
the position of zeros in the structure of mass matrix. Texture six zero mass matrices are 
obtained when both Mi and M v d are texture three zero type, implying that there will be 
400 combinations of texture six zero lepton mass matrices. As mentioned earlier, in the 
present work, we have considered only those mass matrices which lead to non zero and 
distinct mass eigenvalues, therefore imposing the trace and determinant condition on the 
mass matrix, i.e. Trace M^ v b ^ and Det M^ V £> ^ 0, we are left with 12 patterns, clas- 
sified into 2 distinct classes, depending upon the diagonalization equations these satisfy, 
as given in Table ((!]). Details of the diagonalization equations for these 12 mass matrices 
can be looked up in our earlier work [17]. Matrices Mi and M„d each can correspond to 
any of the 12 patterns, therefore yielding 144 possible combinations of texture six zero 
lepton mass matrices which in principle can yield neutrino mixing matrix. These 144 
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Table 1: Table showing various patterns of texture three zero mass matrices categorized 
into two classes, class I and II. 



combinations further form 4 different categories e.g.; 
Category 1: Mi from class I and M vD from class I. 
Category 2: Mi from class II and M vD from class II. 
Category 3: Mi from class I and M v d from class II. 
Category 4: M% from class II and M uD from class I. 

Each category corresponds to 36 possibilities which need exhaustive analysis. As men- 
tioned earlier, we have considered the cases of normally hierarchical, inversely hierarchical 
and degenerate Majorana as well as Dirac neutrinos. Thus, for each of 144 combinations, 
there are 6 cases which have to analyzed leading to a total of 864 possibilities for texture 
six zero lepton mass matrices. 
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5 Majorana Neutrinos 



To begin with, we have analyzed 144 combinations of texture six zero lepton mass ma- 
trices corresponding to normal hierarchy of Majorana neutrinos, by confronting their 
corresponding mixing matrix against latest neutrino oscillation data given in section ([3]). 
It may be mentioned again that for Majorana neutrinos, the hermitian texture three zero 
structure is imposed on M uD as given in equation (J2J). 

In literature, the details of Fritzsch like texture six zero lepton mass matrices [T2] - 
[TB] have been presented by several authors, however, for non-Fritzsch like mass matrices 
similar details have not been presented. In this context, to detail the methodology, we 
consider a typical combination corresponding to category 4, where both Mi and M v d are 
non Fritzsch type, for example, 



a 











A 


A 





a; 









M l = A Ax , M vD = C v B u . (19) 



The diagonalization transformations for these matrices can be easily obtained in terms 
of neutrino masses m Ul , m V2 and m Vi for both normal and inverted hierarchy and are given 
in the Appendix A. The phase matrix Q\P v r> for this particular combination is given as 
diag(e 4flil , e 1 ^ 2 , 1). The scanned ranges of lightest neutrino mass m Ul and phases 0i and 
02 have been given in section ([3]) and AmJ 2 , Am^ and the mixing angles have been 
constrained as given by equations ( 1T7|) and (fl8|) . Using equation ( !T6|) . the PMNS matrix, 
U, obtained for this particular combination is given as 

0.52-0.61 

U = | 0.39 - 0.49 0.40 - 0.53 0.72 - 0.79 | . (20) 

0.67-0.72 

On comparing equation ( |20l) with the ranges of mixing matrix element given by Garcia 
et al. [22] at 3cr C.L., for example, 

0.79-0.86 

U= | 0.25 - 0.53 0.47- 0.73 0.56 -0.79 |, (21) 
0.21 - 0.51 

one finds that the two have good overlap. To graphically show the viability of this 
particular combination of texture six zero mass matrices, in Figure (1) we have plotted 
the lightest neutrino mass against the mixing angle 9i2, #13 and 623 by giving full variation 
to input parameters. The dashed lines depict the limits obtained by assuming normal 
hierarchy and the solid horizontal lines show the experimental 3cr limits of the plotted 
mixing angles. A general look at the Figure (1), shows that mixing angles are well within 
their experimental ranges for a common neutrino mass range for normal hierarchy. 

In a similar manner, one can check the viability of the above set of mass matrices 
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for inversely hierarchical Majorana neutrinos. In the same Figure (1), using dot-dashed 
lines, we have plotted the limits of mixing angles obtained by assuming inverted hierarchy 
against the lightest neutrino mass. It is immediately clear from Figures (la) and (lc) that 
the obtained ranges of 6' 12 and 623 are far from their experimental limits, thus ruling out 
the inverted hierarchy for this particular combination. 

Coming to the degenerate scenarios of Majorana neutrinos characterized by either 
m ui < raj,.-, ~ m U3 ~ 0.1 eV or m vz ~ m Ul < m U2 ~ 0.1 eV corresponding to normal 
hierarchy and inverted hierarchy respectively, one can easily infer from Figures (la) and 
(lc) that degenerate scenario is also ruled out. This can be understood by noting that 
around O.leV the limits obtained by assuming normal hierarchy and inverted hierarchy 
have no overlap with experimental limits of 6*12 and 6*23. 

The viability of rest of the 143 combinations can similarly be checked for normal hier- 
archy, inverted hierarchy and for degenerate Majorana neutrinos. For normal hierarchy, 
we find that in category 1, there are 12 viable combinations. Numerical results corre- 
sponding to one such parallel combination I a I a are given in the first row of Table (2). The 
spectrum of neutrino masses shows that neutrinos are following normal hierarchy. Fur- 
ther, #12 and #13 are spanning their full experimental range, however the obtained range 
of 623 is just below its maximal value. One finds that #23 is sensitive to variations in the 
mass squared differences, however it is not possible to obtain a higher value for 623 even 
when the ranges of Amf 2 and Am^ are extended further. 

Apart from mixing angles, we have also calculated the Jarlskog's rephasing invariant 
in the leptonic sector J;, Dirac like CP violating phase 5i and effective neutrino mass (m ee ) 
related to neutrinoless double beta decay /3/3oi/- The parameter J; has been calculated by 
using the expression 

Ji = lm{U 2Z U* z U 22 UU, (22) 

where U23, U33, U22 and U32 are the elements of mixing matrix U given in equation f )16p . 
The Dirac like CP violating phase 5i can be determined from 

Jl = Sl2S23Sl3Cl2C23C? 3 sin5/, (23) 

where su, S13 and S23 are sines of mixing angles 612, 613 and 623. 

The effective Majorana mass to be measured in f3j3o v decay experiment, is given as 

(m ee ) = m vl Xj\ x + m u2 Ul 2 + m u3 U^ 3 . (24) 

The obtained ranges of Ji, 5i and (m ee ) are in agreement with the ranges obtained in 
other such analysis |16[ [22] . It may be mentioned that the other parallel combinations 
e.g. LjLj, I C I C , 1^1^, I e I e and 1/1/ are isomeric to I a I a , i.e., these combinations are although 
structurally different from each other, their predictions for lepton masses, mixing matrix 
elements and for other phenomenological quantities mentioned above are exactly same 
and hence are not given in the table. Similarly, category 1 also has 6 non parallel viable 
combinations such as I a Ib, isomeric to combinations Id a , IT/, I/I c , Idl e , I e Id, given in 
the second row of the Table (2). It is clear from the table that the results for the non 
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Table 2: Calculated ranges of neutrino masses, mixing angles, (m ee ), Ji and 5\ for viable 
combinations of Mi and M V £> for normally hierarchical Majorana neutrinos. 



parallel combinations are similar to the parallel combinations except that the obtained 
range of 6*23 is above its maximal value. Further, categories 2 and 3 do not have any 
viable combination, because of only two generation mixing in the neutrino mass matrix. 
Lastly, category 4 does have four viable combinations such as II /h isomeric to II^ and 
Ildh isomeric to II/I a , given respectively in the third and fourth row of the Table (2). 
The above two sets can be distinguished again on the basis of 623 , as combination Il/Ib 
gives 6*23 below the maximal value while Il^Ife gives #23 above its maximal value. All the 
four combinations of category 4 lead to a very constrained range of 6*13 e.g. 7° — 9°. 
It may be noted that although most of the phenomenological implications of the above 
mentioned 16 texture six zero lepton mass matrices are similar, however still these can be 
experimentally distinguished with more precise measurements of #23 and #13. The above 
mentioned texture combinations are found to be compatible with the current data even 
when the inputs are at their 2a C.L.. 

In the case of inverted hierarchy and degenerate scenario of Majorana neutrinos, a 
similar analysis has been carried out for all the 144 combinations to check for their com- 
patibility with the latest data. Interestingly, we find that inverted hierarchy as well as 
degenerate Majorana neutrinos are completely ruled out for hermitian texture six zero 
lepton mass matrices. 
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Figure 1: Plots showing variation of mixing angles # 12 , 6> 13 and 6> 2 3 with lightest neutrino 
mass for Ilalb combination given in equation ffl9|) for Majorana neutrinos. The dashed 
lines and the dot-dashed lines depict the limits obtained assuming normal and inverted 
hierarchy respectively, the solid horizontal lines show the experimental 3a limits. 

6 Dirac Neutrinos 

Coming to the case of Dirac neutrinos, we have again analyzed 144 combinations for 
normal hierarchy, inverted hierarchy and degenerate neutrinos. To illustrate the compar- 
ison with the Majorana neutrino case, we pick up the same non-Fritzsch like combination 
given in equation (fl9j) . to check its compatibility with the latest neutrino mixing data. 
The diagonalization transformations for these matrices can easily be obtained in terms of 
neutrino masses m Ul , m U2 and m U3 for both normal and inverted hierarchy and are given 
in the Appendix A. The phase matrix QiP v d and the scanned ranges of lightest neutrino 
mass and phases </>i and 2 have already been mentioned above. 

To check the compatibility of this particular combination for normal hierarchy of Dirac 
neutrinos, in Figure (T5]) we have plotted the allowed parameter space for 8i 2 and 9 2 3 in m Ul - 
9i3 plane, represented respectively by dots and crosses. One can immediately find from 
the figure that there is no common parameter space available to 6' 12 and 6> 23 , concluding 
that this particular combination is not viable for Dirac neutrinos. This result remains 
unaffected even if the input parameter ranges are extended further. Also one can easily 
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check that mixing angle #13 is coming out to be very small. Thus, normally hierarchical 
Dirac neutrinos are ruled out for the texture combination given in equation ffl9|) . Similarly, 
one can also show that degenerate and inversely hierarchical Dirac neutrinos for this 
particular combination are also ruled out. 

The combinations which are viable for Majorana neutrinos, are not viable for Dirac 
neutrinos primarily because of mixing angle #23- For example, for parallel combinations 
given in the first row as well as the combinations given in the third row of the Table (2), 
823 for Dirac neutrinos comes out to be below the experimental limits. Similarly, for non 
parallel combinations given in the second row and the combinations given in the fourth 
row of the Table (2), 623 lies above its experimental limits. 

A similar analysis carried out for the rest of the 128 combinations shows that there 
are no viable texture six zero lepton mass matrices for normally hierarchical, inversely 
hierarchical as well as for degenerate Dirac neutrinos, thus ruling out Dirac neutrinos 
completely for texture six zero mass matrices. 

For the sake of completion, we have also analyzed the cases corresponding to charged 
leptons being in the flavor basis for Dirac as well as Majorana neutrinos and one finds 
that none of matrices gives result within the experimental ranges. 
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Figure 2: Allowed parameter space of 6> 12 and 6*23 in 'lightest neutrino mass - 6*13 plane' 
for normal hierarchy of Dirac neutrinos for texture six zero combination given in equa- 
tion ([19]) . The dots and crosses represent the allowed parameter space for 812 and 823 
respectively. 
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7 Summary and Conclusions 



To summarize, we have analyzed 144 possibilities corresponding to hermitian texture six 
zero lepton mass matrices to check for their viability against current neutrino oscillation 
data. For each of these mass matrices, various cases have been considered in the analy- 
sis, for example, normal hierarchy, inverted hierarchy and degenerate neutrinos for both 
Majorana as well as Dirac neutrinos. For Majorana neutrinos with normal hierarchy, out 
of 144, only 16 combinations are compatible with current neutrino oscillation data at 3a 
C.L.. The above mentioned texture combinations are found to be compatible with current 
data even at 2a level. The 16 viable combinations can be grouped in four sets such that 
the matrices placed in each set are isomeric to each other with regards to their predictions 
for lepton masses and flavor mixing parameters. It is important to note that these four 
sets can be experimentally distinguished from each other with more precise measurements 
of 6*23 and #13. The ranges of neutrino masses, the PMNS matrix, Jarlskog's rephasing 
parameter J 1; Dirac like CP violating phase S\ and effective neutrino mass (m cc ), calcu- 
lated for viable combinations are in agreement with the ranges obtained in other such 
analysis. We find that inverted hierarchy and degenerate neutrinos are completely ruled 
out for texture six zero Majorana neutrinos. Interestingly, for Dirac neutrinos none of the 
144 combinations is viable for normal hierarchy, inverted hierarchy as well as degenerate 
neutrinos, thus ruling out Dirac neutrinos completely for hermitian texture six zero lepton 
mass matrices. 

Acknowledgments 

N.M. would like to thank UGC, Govt, of India, for financial support and the Chairman, 
Department of Physics for providing facilities to work in the department. M.R. and P.S. 
Gill would like to thank the Director, UIET and the Principal, SGGS College, respectively, 
for providing facilities to work. 

References 

[1] R. Davis, Prog. Part. Nucl. Phys. 32, 13 (1994); B. T. Cleveland et al, Astrophys. 
J. 496, 505 (1998); J. N. Abdurashitov et al, SAGE Collaboration, Phys. Rev. C 
60, 055801 (1999); C. M. Cattadori, GNO Collaboration, Nucl. Phys. Proc. Suppl. 
110, 311 (2002); W. Hampel et al, GALLEX Collaboration, Phys. Lett. B 447, 127 
(1999); S. Fukuda et al, Super-Kamiokande Collaboration, Phys. Lett. B 539, 179 
(2002); Q. R. Ahmad et al, SNO Collaboration, Phys. Rev. Lett. 89, 011301 (2002); 
S. N. Ahmad et al, Phys. Rev. Lett. 92, 181301 (2004). 

[2] Y. Fukuda et al, SuperKamiokande Collaboration, Phys. Rev. Lett. 81, 1562 (1998); 
A. Surdo, MACRO Collaboration, Nucl. Phys. Proc. Suppl. 110, 342 (2002); M. 
Sanchez, Soudan 2 Collaboration, Phys. Rev. D 68, 113004 (2003). 

[3] K. Eguchi et al, KamLAND Collaboration, Phys. Rev. Lett. 90, 021802 (2003); 
Phys. Rev. Lett. 94, 081801 (2005). 



12 



[4] M. H. Aim et ai, K2K Collaboration, Phys. Rev. Lett. 90, 041801 (2003). 

[5] H. Fritzsch, Phys. Lett. B 70, 436 (1977); B 73, 317 (1978); S. Weinberg, in Trans- 
actions of the New York Academy of Sciences 38, 185 (1977); F. Wilczek and A. Zee, 
Phys. Lett. B 70, 418 (1977). 

[6] P. Minkowski, Phys. Lett. B 67 421 (1977); T. Yanagida, in Proceedings of the 
Workshop on Unified Theory and the Baryon Number of the Universe, edited by O. 
Sawada and A. Sugamoto (KEK, Tsukuba, 1979), p. 95; M. Gell-Mann, P. Ramond, 
and R. Slansky, in Supergravity, edited by F. van Nieuwenhuizen and D. Freedman 
(North Holland, Amsterdam, 1979), p. 315; S.L. Glashow, in Quarks and Leptons, 
edited by M. Levy et al. (Plenum, New York, 1980), p. 707; R. N. Mohapatra and 
G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[7] S. Weinberg, Phys. Rev. D 5, 1962 (1972); A. Zee, Phys. Lett. B 93, 389 (1980). 

[8] CD. Froggatt and H.B. Nielsen, Nucl. Phys. B 147, 277 (1979); H. Harari, H. Haut, 
and J. Weyers, Phys. Lett. B 78, 459 (1978). 

[9] K. R. Dienes, E. Dudas, and T. Gherghetta, Nucl. Phys. B 557, 25 (1999); N. 
ArkaniHamed, S. Dimopoulos, G. Dvali, J. March- Russell, Phys. Rev. D 65, 024032 
(2002)( p^ph/9811448D ; G. Dvali and A.Yu. Smirnov, Nucl. Phys. B 563, 63 (1999). 

[10] H. Fritzsch, Z. Z. Xing, Prog. Part. Nucl. Phys. 45, 1 (2000), and references therein. 

[11] P. S. Gill, M. Gupta, J. Phys. G: Nuc. Part. Phys. 21, (1995) 1; P. S. Gill, M. Gupta, 
Phys. Rev. D 57, 3971 (1998) ; M. Randhawa, V. Bhatnagar, P. S. Gill, M. Gupta, 
Phys. Rev. D 60, 051301 (1999), M. Randhawa, M. Gupta, Phys. Rev. D 63, 097301 

(2001) . 

[12] H. Fritzsch, Z. Z. Xing, Int. J. Mod. Phys. A 19, 1 (2004) and references therein. 

[13] K. Kang, S. K. Kang, Phys. Rev. D 56, 1511 (1997); H. Nishiura, K. Matsuda, T. 
Fukuyama, Phys. Rev. D 60, 013006 (1999); K. Matsuda, T. Fukuyama, H. Nishiura, 
Phys. Rev. D 61, 053001 (2000) ; K. Kang, S. K. Kang, C. S. Kim, S. M. Kim, Mod. 
Phys. Lett. A 16, 2169 (2001); C. Giunti, M. Tanimoto, Phys. Rev. D 66, 113006 

(2002) ; M. Frigerio, A. Yu. Smirnov, Nucl. Phys. B 640, 233 (2002); P. F. Harrison, 
W. G. Scott, Phys. Lett. B 547, 219 (2002) ; E. Ma, Mod. Phys. Lett. A 17, 2361 
(2002); E. Ma, Phys. Rev. D 66, 117301 (2002); G. Bhattacharyya, A. Raychaudhuri, 
A. Sil, Phys. Rev. D 67, 073004 (2003); P. F. Harrison, W. G. Scott, Phys. Lett. 
B 594, 324 (2004); M. Bando, S. Kaneko, M. Obara, M. Tanimoto, Phys. Lett. B 
580, 229 (2004); O. L. G. Peres, A. Yu. Smirnov, Nucl. Phys. B 680, 479 (2004); 
C. H. Albright, Phys. Lett. B 599, 285 (2004); J. Ferrandis, S. Pakvasa, Phys. Lett. 
B 603, 184 (2004); S. T. Petcov, W. Rodejohann, Phys. Rev. D 71, 073002 (2005); 
C. Hagedorn, W. Rodejohann, JHEP 0507, 034 (2005); S. S. Masood, S. Nasri, J. 
Schechter, Phys. Rev. D 71, 093005 (2005); R. Derrmsek, S. Raby, Phys. Lett. B 
622, 327 (2005); F. Plentinger, W. Rodejohann, Phys. Lett. B 625, 264 (2005). 



13 



[14] P. H. Frampton, S. L. Glashow, D. Marfatia, Phys. Lett. B 536, 79 (2002); Z. Z. 
Xing, Phys. Lett. B 530, 159 (2002); Z. Z. Xing, Phys. Lett. B 539, 85 (2002); A. 
Kageyama, S. Kaneko, N. Simoyama, M. Tanimoto, Phys. Lett. B 538, 96 (2002); 
W. L. Guo, Z. Z. Xing, |hep-ph / 02 1 1 3 1 5 [ M. Frigerio, A. Yu. Smirnov, Phys. Rev. 
D 67, 013007 (2003); K. S. Babu, E. Ma, J. W. F. Valle, Phys. Lett. B 552, 207 
(2003); W. L. Guo, Z. Z. Xing, Phys. Rev. D 67, 053002 (2003); S. Kaneko, M. 
Tanimoto, Phys. Lett. B 551, 127 (2003); B. R. Desai, D. P. Roy, A. R. Vaucher, 
Mod. Phys. Lett. A 18, 1355 (2003); K. Hasegawa, C. S. Lim, K. Ogure, Phys. Rev. 
D 68, 053006 (2003); M. Honda, S. Kaneko, M. Tanimoto, JHEP 0309, 028 (2003); 
Z. Z. Xing, H. Zhang, Phys. Lett. B 569, 30 (2003); C. Hagedorn, W. Rodejohann, 
JHEP 0507, 034 (2005); E. I. Lashin, N. Chamaon, Phys. Rev. D 78, 073002 (2008); 
E. I. Lashin, N. Chamaon, Phys. Rev. D 80, 093004 (2009). 

[15] S. Zhou, Z. Z. Xing, Eur. Phys. J. C 38, 495 (2005). 

[16] M. Randhawa, G. Ahuja, M. Gupta, Phys. Rev. D 65, 093016 (2002); Z. Z. Xing, 
H. Zhang, Phy. Lett B 569, 30 (2003); Z. Z. Xing, S. Zhou, Phy. Lett B 593, 
156 (2004); M. Randhawa, G. Ahuja, M. Gupta, Phys. Lett. B 643, 175 (2006); G. 
Ahuja, S.Kumar, M. Randhawa, M. Gupta, S. Dev, Phys. Rev. D 76, 013006 (2007); 
S. Choubey, W. Rodejohann, P. Roy, Nucl. Phys. B 808, 272 (2009), Erratum-ibid. 
818:136, 2009; H. Fritzsch, X. X. Xing, Phys. Lett B 682, 220 (2009), S. Goswami, 
A. Watanabe, Phys. Rev. D 79, 033004 (2009); G. Ahuja, M. Gupta, M. Randhawa, 
R. Verma, Phys. Rev. D 79, 093006 (2009). 

[17] N. Mahajan, R. Verma, M. Gupta, Int. J. Mod. Phys. A 25, 1 (2010). 

[18] A. Strumia, F. Vissani, |hep-ph/0606054l 

[19] Michael A. Schmidt, Alexei Yu. Smirnov, Phys. Rev. D 74, 113003 (2006). 

[20] B. Pontecorvo, Zh. Eksp. Teor. Fiz. (JETP) 33, 549 (1957); ibid 53, (1967) 1717; Z. 
Maki, M. Nakagawa, S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 

[21] M. C. Gonzalez-Garcia, M. Maltoni, J. Salvado, JHEP 04, 056 (2010). 

[22] M. C. Gonzalez-Garcia, M. Maltoni, Phys. Rept. 460, 1 (2008). 



14 



Appendix A 

The diagonalization transformations for the real texture 3 zero mass matrices given in 
equation ffl9l) are as follows: 



1. The diagonalizing matrix Oi for the real 3 zero mass matrix Mi, 




(25) 



is given as 



/ 1 



O, 



\ 



(26) 



m M +m T 

2. The diagonalizing matrix O v d for the real 3 zero neutrino mass matrix M u o, 



M, 



uD 









A v 





c v 


B v 


K 


Bl 






(27) 



for normal hierarchy of neutrinos is given as 



m 2 m 3 (m 3 —m2) 



\ 



(mj -m 2 +m 3 )(mi+m 2 )(m 3 -mi) 

mi (mi + m 3 ) (m 2 —mi) 

( m i — m 2 + m 3 ) ( m 1 + m 2 ) ( m 3 — mi) 
("T-3 -m 2 ) 



mim 3 (mi+m 3 ) 



(mi — m 2 + m 3 ) (mi +m 2 ) (m3+m 2 ) 

m 2 (m 2 -mi)(m3 — m 2 ) 

(ni 1 m 2 ■ ; . ■■■■■ j + m 2 )(m 2 + m 3 ) 
m 2 (mi + m 3 ) 



("il - 



m 2 mi(m 2 — mi) 



V (mi+m 2 )(m 3 — mi) y ( m l +m 2 ) (m 3 +m 2 ) 

Similarly, for inverted hierarchy of neutrinos, O v u is given as 



m-2 /7i,;;(m.;:;+m2 ) 



o v 



\ 



( — mi +m 2 +m 3 )(mi+m 2 )(m 3 +mi) 
/ mi (mi — m 3 )(m 2 — mi) 

/ ( — mi +m 2 +m 3 )(m.i+m 2 )( m 3 — mi) 
/ mi(m 3 +m 2 )~ 

V ( 



mim 3 (mi -m 3 ) 



(mi+m 2 )(m 3 +mi) 



/ ( — mi4 fn j m 3 Kmi ~-ffl 2 )(m3+m 2 ) 
/ m 2 (m2 — m 1 )(m3 + TTi 2 ) 

( — m 1 + m 2 +m3)(m 1 +m 2 )(m 2 — m 3 ) 
m 2 (mi— m 3 ) 



(mi +m 2 )(m 2 — 7713) 



m 2 +m 3 )(m3+m 2 )(m3-m 1 ) 

m3(m 1 +m 3 )(m3-m 2 ) 
1"! -i Ti i |m 3 )(m i |m2)(??i i -mi) 

V (">3-n>l)(n>3+ m j) 



/ m 2 mi(mi— m 2 ) 

( — m 1 + m 2 + m 3 ) ( n 1 ; ; - 7 n 1 ) ) ; 3 — m 2 ) 
/ m 3 (m 2 +m3)(m 3 -mi) 
/ (-mi +m 2 + m3)(m3-m 2 )( m 3 +mi) 
[ m 3 (mi-m 2 ) 
V (m 3 -m 2 )(m3+mi) 



where mi 



m,. 



m 2 



m,. 



and m 3 



for Dirac neutrinos and m\ = y/m Ul m^, 



m 2 = y , rn V 2 m R an d ^3 = y/m U3 mji for Majorana neutrinos. 
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